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ABSTRACT 

Polycrystalline  alloys  of  iron  in  vanadium  up  to  15.32 
atomic  percent  were  prepared  principally  by  arc  melting.   In 
order  to  obtain  the  elastic  constants  and  investigate  the 
possibility  of  the  existence  of  an  elastically  isotropic 
alloy  ultrasonic  measurements  were  made.   These  measurements 
determined  the  longitudinal  and  transverse  wave  velocities 
in  the  10-30  MHz  range  at  room  temperature.   The  resulting 
longitudinal  modulus  increased  with  alloying  whereas  the 
shear  modulus  remained  roughly  constant  although  the  data 
are  somewhat  scattered.   These  data  are  interpreted  to  show 
that  the  shear  constants  vary  through  alloying  in  a  manner 
consistent  with  the  results  of  other  experiments  and  that 
the  results  are  consistent  with  a  s1  d      electronic 
configuration  in  the  transition  metals.   The  existence  of 
an  isotropic  alloy  is  indicated  for  alloys  with  electron  to 
atom  ratios  lower  than  vanadium  rather  than  for  the  alloys 
measured. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 5 

II.  EXPERIMENTAL  DETAILS 8 

A.  SAMPLE  PREPARATION 8 

B.  MEASUREMENT  OF  THE  SOUND  VELOCITY 9 

III.  EXPERIMENTAL  RESULTS 15 

IV.  DISCUSSION 18 

V.  CONCLUSIONS 2  8 

TABLES  I  TO  III 30 

FIGURES  1  TO  6 33 

APPENDIX  I 39 

APPENDIX  II 4  8 

REFERENCES 56 

INITIAL  DISTRIBUTION  LIST 58 

FORM  DD  1473 59 


ACKNOWLEDGEMENTS 

The  author  wishes  to  acknowledge  the  Bureau  of  Mines 
for  providing  most  of  the  vanadium.   Also  Dr.  Martin 
Glicksman  and  his  group  at  the  Naval  Research  Laboratory 
are  thanked  for  their  assistance  in  preparation  of  the 
alloys.   Finally  much  gratitude  is  owed  to  John  R. 
Neighbours  and  Eugene  C.  Crittenden,  Jr.,  Professors  of 
Physics,  Naval  Postgraduate  School  whose  guidance  made 
this  work  possible. 


I .   INTRODUCTION 

The  decision  to  investigate  the  elastic  properties  of  a 
vanadium- iron  alloy  was  the  result  of  a  desire  to  have  an 
elastically  isotropic  material  such  as  tungsten,  but  without 
tungsten's  high  melting  point,  for  use  in  a  thin  film 
experiment.   It  was  felt  that  the  isotropic  property  of  such 
a  material  would  also  be  advantageous  in  macroscopic  poly- 
crystalline  form.   Isotropic  grains  would  tend  to  deform 
uniformly  in  the  presence  of  an  overall  stress  irregardless 
of  their  orientations.   This  effect  would  reduce  one  source 
of  motion  at  the  grain  boundaries  resulting  in  a  corresponding 
improved  resistence  to  cyclic  fatigue. 

It  was  known  that  the  ratio  of  the  elastic  shear  moduli, 
A  =  C/C ' ,  called  the  Zener  anisotropy  is  less  than  unity  for 
some  of  the  bcc  transition  elements  and  greater  than  unity 
for  others.   Using  a  simple  Vegards ' s  Law  model  for  the 
variation  of  A  with  alloying,  various  phase  diagrams  were 
investigated  to  find  a  pair  of  constituents  where  the  re- 
quired relative  concentrations  for  isotropy  fell  in  a  region 
of  solid  solubility  of  the  system.   The  vanadium-iron  system 
was  found  to  have  a  solid  solution  bcc  phase  up  to  approxi- 
mately 40  at.%  iron  in  vanadium.   On  the  basis  of  Vegard's 
Law  the  concentration  of  iron  for  isotropy  was  estimated  to 
be  approximately  12  at.%.   Another  influencing  factor  in  the 
choice  of  the  V-Fe  system  was  the  recent  availability  of 
very  high  purity  vanadium. 


The  free  atom  electron  configuration  of  vanadium  is 
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4s  3d   and  of  iron,  4s  3d  .   Thus  the  addition  of  iron  to 
vanadium  increases  the  electron  to  atom  ratio  of  the  alloy. 
According  to  Fuchs  [Ref.  1]  the  resistance  to  shear  in  metals 
is  a  result  of  two  interactions:  (a)  the  Coulomb  attraction 
between  positive  ions  and  negative  conduction  electrons,  and 
(b)  the  exchange  interaction  between  ions,  thought  to  be 
significant  only  between  nearest  and  second  nearest  neighbors 
This  second  interaction  depends  in  a  sensitive  way  -  on  the 
electron  configuration  and,  in  transition  metal  alloys  in 
particular,  the  d  shell  electrons  are  thought  to  contribute 
significantly. 

Very  little  work  has  been  done  on  the  elastic  constants 
of  binary  alloys  where  both  constituents  are  transition 
metals.   In  addition  to  making  an  isotropic  alloy  some 
knowledge  about  the  interaction  between  atoms  in  such  an 
alloy  would  be  expected  to  be  gained  from  the  analysis  of  a 
series  of  compositions.   Therefore,  seven  polycrystalline 
compositions  of  the  vanadium-iron  system  were  produced  with 
iron  content  between  3.25  and  15.32   at.%.   The  longitudinal 
and  transverse  elastic  constants  were  determined,  by  a  pulse 
echo  method,  for  each  composition.   Then,  using  an  empirical 
relation  for  the  shear  constant  C'=i/2 (C   -C   ),  the  single 
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crystal  elastic  constants  were  calculated  and  discussed  in 
light  of  the  variations  in  electron  to  atom  ratio  and  the 
d-band  contributions. 

The  classical  theory  of  elasticity  is  discussed  in 
Appendix  I  and  the  Fuch's  theory  of  determination  of  the 


elastic  constants  from  fundamental  considerations  is  discussed 
in  Appendix  II. 


II.   EXPERIMENTAL  DETAILS 

A.   SAMPLE  PREPARATION 

The  preparation  of  the  seven  alloys  used  in  this  study 
was  done  at  the  Naval  Research  Laboratory,  Washington,  D.C. 
The  history  of  each  sample  measured  is  shown  in  Table  I. 
Commercial  analysis  of  the  composition  and  impurities  of 
each  sample  is  shown  in  Table  II. 

Two  methods  of  producing  the  alloy  were  used.   The 
constituents  were  melted  by  an  electron  beam  in  vacuum  or 
arc-melting  under  an  ultra-pure  helium  atmosphere.   Electron 
beam  melting  was  discontinued  after  the  second  sample  when 
it  was  determined  that  the  refining  effect  of  this  method 
reduced  the  iron  content  of  the  alloy  by  unpredictable 
amounts.   Ingots  from  both  methods  were  hot  forged  and 
swaged  into  rods.   A  typical  rod  was  3/8  to  1/2  inch  in 
diameter  and  12  to  15  inches  long. 

These  rods  were  used  in  attempts  to  grow  alloy  single 
crystals.   The  attempted  growth  methods  were  liquid  float 
zoning  and  strain  annealing.   All  of  the  samples  exhibited 
grain  growth:  zone  melting  produced  grains  of  .25  to  .75 
mm  diameter,  while  strain  anneal  methods  produced  grains 
of  1.25  to  3.75  mm  diameter.   Time  and  the  availability  of 
material  precluded  a  more  thorough  exploration  of  the  pos- 
sibilities of  producing  single  crystals  and  thus  poly- 
crystalline  rods  of  each  composition  were  selected  for 
measurement. 
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Fourteen  specimens  were  made  from  the  seven  alloys. 
Seven  were  2  cm  long  and  .6  35  cm  in  diameter.   These  were 
used  for  measuring  the  longitudinal  velocity.   The  remain- 
der, cut  adjacent  to  the  first,  were  1  cm  long  and  were 
used  for  measuring  the  transverse  velocity.   To  provide 
the  best  possible  diameter  to  wavelength  ratio  in  the 
transverse  samples  the  largest  available  diameter  was  pre- 
served in  each  case.   Sample  diameters  varied  from  .6  35  cm 
to  1.10  cm.   The  ends  of  all  of  the  cylindrical  specimens 
were  machined  flat  and  parallel  and  the  surface  polished. 
The  surfaces  were  found  to  be  parallel  to  an  average  of 
1.5  milliradians . 

B.   MEASUREMENT  OF  THE  SOUND  VELOCITY 

The  velocity  of  sound  in  the  seven  alloys  was  measured 
by  the  pulse  echo  technique.   A  sound  pulse  was  transmitted 
into  the  material  by  a  piezo-electric  quartz  transducer 
driven  by  a  pulse  generator  operating  at  the  transducer 
fundamental  frequency  or  an  odd  harmonic.   The  initial 
pulse  and  each  succeeding  echo  can  be  received  by  a  second 
transducer.   This  second  transducer  is  placed  at  the  oppo- 
site end  of  the  sample.   For  the  transverse  measurement  the 
two  transducers  must  be  aligned  so  that  their  particle 
motions  are  parallel.   The  two  transducer  method  eliminated 
saturation  of  the  receiver  amplifier  by  the  driving  pulse. 
The  arrival  time  of  each  echo  is  usually  determined  by 
using  an  oscilloscope  trace  with  a  delay  time  multiplier. 


The  transducers  used  were  10MHz  X  cut  quartz  for  the 
longitudinal  mode  measurements  and  AC  cuts  of  the  same 
frequency  for  the  transverse  mode.   The  transducers  were 
bonded  to  the  samples  with  Salol.   A  spring  loaded  jig  was 
used  so  that  both  transducers  could  be  bonded  simultaneously. 
This  improved  the  uniformity  of  the  bonds  and  eliminated 
the  problem  of  melting  the  first  bond  while  the  sample  was 
being  heated  for  the  second.   During  the  measurements  the 
sample  was  suspended  at  its  center  by  a  three  point  circular 
holder  which  was  concentric  with  two  shielded  contacts  spring 
loaded  to  the  transducers.   Series  coils  were  employed  to 
provide  a  low  impedance  match  of  the  capacitive  transducers 
to  the  pulse  generator  and  receiver.   All  measurements  were 
performed  at  ambient  room  temperature  (29  4  K) . 

The  alloys  in  this  experiment  were,  in  most  cases,  highly 
attenuating.   The  transverse  mode  was  generally  more  at- 
tenuated than  the  longitudinal  mode  and  often  only  two  or 
three  echos  could  be  discerned  above  the  system  noise. 
Measurements  of  the  total  attenuation  were  made  in  the 
longitudinal  case.   The  total  attenuation  includes  all  loses, 
which  are  loses  at  the  bonds  and  transducers  as  well  as  in 
the  sample  themselves.   These  loses  varied  from  a  low  of 
3db  per  cm  for  the  best  sample  to  approximately  5  db  per  cm 
for  the  worst  sample.   For  this  reason  the  unrectified  pulse 
waveform,  after  amplification,  was  further  processed  by  a 
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Princeton  Applied  Research  (PAR)  Model  160  Boxcar 
Integrator. 

The  boxcar  integrator  has  three  principal  advantages  for 
the  pulse  echo  experiment.   It  provides  a  large  improvement 
in  signal  to  noise  ratio.   It  faithfully  reproduces  the  small 
amplitude  portions  of  the  waveform  so  that  the  arrival  time 
of  the  pulse  can  be  accurately  determined.   In  conjunction 
with  a  strip  chart  recorder  it  provides  a  permanent  working 
data  record.   The  boxcar  integrator  operates  by  scanning  an 
aperture  in  time  along  the  waveform  and  averaging  the  value 
of  the  signal  at  each  point  over  many  repititions  of  the 
waveform.   The  result  is  that  the  random  noise  fluctuations 
on  the  signal  will  average  to  zero,  and  if  enough  repetitions 
of  the  signal  are  integrated,  the  output  will  be  the  true 
value.   Typically  50  usee  of  the  echo  train  was  scanned 
using  a  10  nsec  aperture.   Such  an  aperture  can  resolve 
frequencies  up  to  about  30  MHz.   The  scanning  time  to  make 
each  pass  along  the  train  was  50  minutes.   A  pulse  repitition 
frequency  of  10  KHz  was  used  so  that  each  point  on  the  wave- 
form was  sampled  6x10 3  times.   The  output  to  the  strip  re- 
corder is  then  a  picture  of  the  echo  train  produced  at  one 
microsecond  per  minute.   By  diverting  the  boxcar  input  to 
the  time  mark  generator  that  was  triggering  the  experiment 
a  trace  of  convenient  time  marks  could  be  placed  in  the 
"dead  space"  between  echos.   These  time  marks  could  not  be 
used  to  measure  the  arrival  time  of  echos,  however,  since 
any  variation  in  the  linearity  of  the  travel  of  the  aperture 
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along  the  signal  (ie  along  the  timebase)  or  any  variation 
of  the  speed  of  the  strip  chart  recorder  would  have  the 
effect  of  expanding  or  contracting  the  trace.   Instead  the 
voltage  corresponding  to  the  position  of  the  aperture  along 
the  timebase  was  used  to  drive  a  second  pen  which  ran  con- 
currently with  the  echo  trace.   (Figure  1)   This  second 
trace  was  used  to  determine  the  echo  arrival  times.   The 
accuracy  of  the  time  measurement  is  improved  by  this  re- 
finement since  if  the  nonlinearities  above  occur  they  will 
be  automatically  compensated.   The  time  marks  on  the  echo 
trace  were  used  to  check  the  calibration  of  the  time  trace. 
Electronic  instability  over  the  relatively  long  times  re- 
quired to  make  a  trace  did  not  present  a  problem.   Varia- 
tions in  the  trigger  rate  of  the  time  mark  generator  or  in 
the  trigger  response  of  the  pulse  generator  and  boxcar 
integrator  that  would  have  resulted  in  noise  in  the  output 
trace  were  not  observed  to  occur.   From  these  traces  the 
arrival  time  of  each  echo  was  measured  and  plotted  against 
the  elapsed  path  length.   The  slope  of  this  plot  is  the 
sound  velocity. 

In  a  homogeneous  polycrystalline  material  the  value  of 
the  transverse  velocity  should  not  depend  on  the  orientation 
of  the  transducer  with  respect  to  the  sample.   Such  is  not 
the  case  in  single  crystal  measurements  where  alignment  of 
the  particle  motion  with  certain  crystallographic  directions 
is  required  to  determine  the  elastic  constants.   The  rather 
large  grain  size  in  these  samples  could  have  the  effect  that 
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the  orientations  of  the  grains  might  not  be  sufficiently 
random  to  result  in  the  true  average  velocity  being  measured 
This  effect  of  a  lack  of  randomness  is  reduced  if  the  aniso- 
tropy  of  the  individual  grains  is  low,  and  vice  versa.   To 
compensate  for  any  lack  of  random  grain  orientation  velocity 
measurements  were  made  on  each  sample  with  the  induced 
particle  motion  oriented  at  0,  60  and  120  degrees  from  an 
arbitrary  mark.   The  velocities  thus  measured  and  several 
values  at  unknown  orientations  were  combined  by  taking  the 
arithmetic  mean  as  the  velocity  for  the  composition.   The 
choice  of  a  difference  in  angle  of  11/ 3  for  the  various 
orientations  was  chosen  as  convenient  and  because  a  cubic 
material  has  II /^  symmetry.   The  variation  of  sound  velocity 
with  orientation  for  each  sample  is  indicated  by  the  bars 
on  the  data  in  Figure  4 .   These  bars  are  drawn  between  the 
maximum  and  minimum  velocities  observed.   The  average  varia- 
tion in  velocity  between  orientations  and  samples  was  2.26%. 
The  repeatability  on  a  given  sample  for  a  given  orientation 
is  within  the  accuracy  of  the  experiment  which  is  discussed 
below. 

To  determine  the  elastic  constants  both  the  sound 
velocity  and  the  density  of  the  sample  must  be  known. 
Pearson  [Ref.  2]  lists  data  by  Hanneman  and  Mariano  which 
includes  the  value  of  the  lattice  parameter  of  the  vanadium- 
iron  system  for  0,  9.5  and  20.5   at.%  iron.   If  a  smooth 
curve  is  drawn  through  these  points  the  lattice  parameter  of 
the  seven  alloy  compositions  used  can  be  found.   The  density 
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determined  by  this  method  is  compared  in  Figure  2  with  the 
density  determined  by  weighing  and  measuring  the  samples. 
Figure  2  shows  density  as  a  function  of  the  atomic  percent- 
age of  iron  and  indicates  the  agreement  between  the  two 
methods  is  within  1%  for  all  compositions.   The  experimental 
points  can  also  be  fitted  to  a  straight  line  with  one  point 
outside  of  1%  and  an  average  deviation  of  0.3%.   Such  a 
straight  line  fit  would  increase  the  density  determination 
of  the  dilute  composition  approximately  0.5%  and  decrease  the 
concentrated  composition  a  like  amount.   It  is  felt  that  the 
precision  of  the  x-ray  determination  of  the  lattice  parameter 
exceeded  that  of  the  experimental  check  on  the  calculated 
density  so  that  the  density  of  each  composition  has  been 
calculated  from  the  lattice  parameter  read  from  a  smooth 
curve  through  the  published  values. 
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III.   EXPERIMENTAL  RESULTS 

The  elastic  constants  of  a  cubic  polycrystalline  material 
are  the  shear  constant  G,  and  the  longitudinal  constant, 
(B  +  —  G) ,  where  B  is  the  bulk  modulus.   Appendix  I  shows 

3 

that  the  longitudinal  and  transverse  sound  velocities  are 
given  by 


V,  =  JB+  -  G 


P 

In  order  to  compare  the  quantities  measured  for  the 
alloys  with  those  of  pure  vanadium  the  Hill  method  [Ref.  3] 
was  used  to  average  over  all  lattice  orientations  the  single 
crystal  elastic  constants  of  vanadium  as  measured  by  Alers 
[Ref.  4].   This  average  is  required  to  obtain  the  poly- 
crystalline value  of  G  and  is  also  discussed  below.   The 
bulk  modulus  in  cubic  materials  is  the  same  for  the  single 
crystal  and  polycrystal  case.   After  a  small  correction 
for  the  difference  between  the  density  used  by  Alers  and 
that  calculated  from  the  lattice  parameter  B  and  G  for  pure 
vanadium  are  known.   For  vanadium  and  each  composition,  the 
atomic  percentage  of  iron,  density,  longitudinal  constant, 
shear  constant  and  bulk  modulus  are  given  in  Table  III. 
Graphs  of  the  measured  elastic  constants  which  also  contain 
the  calculated  values  for  vanadium  are  shown  in  Figure  3 
to  Figure  5. 


15 


The  longitudinal  modulus  (Figure  3)  increases  with  in- 
creasing iron  concentration.   The  experimental  data  is  fitted 
best  by  a  least  squares  parabola.   With  the  exception  of  the 
3.25   at.%  composition,  which  is  approximately  2%  low,  all 
of  the  other  samples  and  the  vanadium  value  fall  within  1% 
of  this  fit.   Above  approximately  5   at.%  iron  the  measured 
values  increase  essentially  linearly.   However,  a  linear 
least  squares  fit  to  the  data  had  a  larger  least  square 
error  and  fell  approximately  3%  low  for  the  value  of  pure 
vanadium. 

The  transverse  modulus  (Figure  4)  is  essentially  con- 
stant over  the  entire  range  of  compositions  investigated. 
These  is  considerably  more  scatter  in  the  data  for  the 
transverse  case  than  in  the  longitudinal  case.   The  best 
fit  straight  line  increases  approximately  1.5%  over  the 
pure  vanadium  value  for  the  15.32   at.%  composition.   All 
of  the  samples  fall  within  5%  of  this  line  with  the  ex- 
ception of  the  7.40   at.%  composition  which  is  approximately 
12%  low.   The  reproducibility  of  this  low  velocity  for 
sample  3  was  demonstrated  on  a  third  specimen  cut  adjacent 
to  the  first  two.   Because  of  the  small  size  of  the  transverse 
constant  with  respect  to  the  longitudinal  constant  the 
percentage  errors  are  larger  although  the  absolute  errors 
are  approximately  the  same. 

The  bulk  modulus  (Figure  5)  is  readily  obtained  from 
the  definition  of  the  longitudinal  constant  (B  +  —  G) . 
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Since  it  is  a  combination  of  both  the  longitudinal  and 
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transverse  constants  it  reflects  both  the  rise  with  in- 
creasing iron  concentration  of  the  longitudinal  constant 
and  the  scatter  present  in  the  transverse  data. 

Based  on  the  estimate  of  the  errors  in  the  measured 

time,  length  and  density  the  accuracy  of  a  single  measure- 

2 
ment  of  pV   is  estimated  to  be  approximately  .75%.   The 

largest  portion  of  this  error  is  in  the  time  determination 

Differences  in  individual  echo  arrival  times  for  different 

runs  on  the  same  sample  that  resulted  in  the  same  velocity 

were  found  to  be  approximately  .75%. 
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IV.   DISCUSSION 

The  two  elastic  constants  measured  in  a  polycrystalline 

experiment  are  the  shear  constant,  G,  and  the  longitudinal 

if 
constant  (B  +  —  G) .   For  this  experiment  therefore  the 

3 

experimental  results  are  the  variation  of  the  shear  modulus, 
G,  and  the  bulk  modulus,  B,  with  increasing  atomic  per- 
centages of  iron. 

Although  the  interpretation  of  these  results  alone  is 
interesting  they  do  not  provide  as  much  information,  that 
would  be  useful  in  interpreting  the  fundamental  contributions 
of  atomic  forces  to  the  elastic  constants,  as  the  single 
crystal  elastic  constants  would.   (See  Appendix  II  for 
the  theory  of  the  fundamental  contributions  to  elastic 
constants)   Moreover,  for  cubic  materials  there  are  three 
non-zero  independent  elastic  constants,  (See  Appendix  I) 
and  it  is  impossible  to  unfold  these  three  quantities 
from  the  two  measured  quantities. 

An  independent  source  of  the  third  constant  was  avail- 
able in  the  work  of  Fisher  and  Dever  on  the  stability  of 
bcc  transition  metals  [Ref .  5] .     They  present  an  empirical 
relation  for  the  C1  elastic  shear  constant  as  a  function 
of  the  electron  to  atom  ratio  of  transition  metal  alloys. 
Here  C1  is  one  of  a  set  of  3  constants  which  are  linear 
combinations  of  the  basic  cubic  constants,  Ci i  C12  and 
Ci*tt  (See  equation  A1.8,f.f.). 


18 


Ci  i  +  2Ci 2 

B  =  


Ci i  -  Ci2  (1) 

C'= 


c  =  c 


»t  t» 


Here  B  is  the  bulk  modulus  and  C  and  C1  are  shear  constants. 
In  terms  of  this  set  of  constants  the  measured  bulk  modulus 
for  the  polycrystalline  material  is  directly  comparable  with 
B=l/3 (Ci i+2Ci 2 ) ,    however,  the  shear  modulus,  G,  is  a  function 
of  both  C  and  C.   The.  fact  that  G  is  a  known  function  of 
C  and  C  is  a  result  of  the  methods  of  determining  poly- 
crystalline aggregate  moduli,  as  averages   of  the  single 
crystal  moduli.   It  is  thus  instructive  to  digress  briefly 
to  discuss  these  averages. 

In  1952  Hill  [Ref.  3]  showed  that  two  classical  approxi- 
mation to  polycrystalline  averages  by  Voigt  in  1928  [Ref.  6] 
and  Reuss  in  1929  [Ref.  7]  were  upper  and  lower  bounds  on 
the  problem  and  suggested  the  arithmetic  mean  as  an  average 
that  would  be  closer  to  measured  values.   Voigt  averaged 
over  all  lattice  orientations  using  the  assumption  that 
the  strain  is  uniform  throughout  a  grain.   Reuss  did  the 
averaging  assuming  that  the  stress  is  uniform  throughout 
a  grain.   Following  Hill,  the  strain  energy  function  of 
an  isotropic  material  in  terms  of  the  principal  stresses 

T.  and  strains  e.  and  the  unknown  bulk  and  shear  modulus, 
i  1 

B  and  G  are 
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3 
2 

w  =  B(ei+e2+e3) 2+  -G    S   (e.-e.)2      (2) 

3   i5,j=i   i   D 

3 
1  3 

w  =  -(T1+T2+T3)2+  -p    £   (T.-T.)2      (3) 

Even  though  the  Voigt  and  Reuss  approximation  are  limits 
they  lead  to  simplified  relations  between  the  polycrystal- 
line  constants  and  the  single  crystal  constants. 


9Bv  =  (Ci 1+C22+C33)  +  2(Ci2+C23+C3i)  (4) 

15Gv  =  (Cii+C22+C33)  -  (Ci2+C23+C3i)  +3  (Cim+Css+Cee)    (5) 

g-  =  (Si i+S22+S33)  +  2(Si2+S23+S3i)  (6) 

r 

^-  =  (Sii+S22+S33)  -4(S12+S23+S3i)  +3  (Sm+Sss+Sgb)     (7) 
Gr 

These  relations  apply  to  all  crystal  classes  and  the  sub- 
scripts v  and  r  deonte  Voigt   and  Reuss  respectively.   The 
S.  .  are  the  elastic  compliances  (see  Appendix  I). 

By  writing  the  C. .  or  S. .  in  (4) through  (7)  in  terms  of 
the  cubic  constants  Cn,  C12  or  Cifi»  and  using  the  cubic 
inverse  relations 

r      _       (n        n      \  _  1 i     (Cn+  2C12)  =  

(S11-S12) 

one  obtains, 

1 

B   =  B   =  -  (C11+2C12) •  (8) 

v    r    3 
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which  is  identical  to  the  single  crystal  experssion  in 
equation  (1) ,  and 

G  =  '(Cn-Cu+SCO  „  5^(Cn-C12)    %  (9) 

v    5  r   4C^-3(Cn-C12) 

Now  converting  to  C  and  C1  using  equation  (1) 

2     3                      10CC1 
G  =  -C,+  -  C  G  =  ±^= .  (10) 

v    5     5  r    4C  +  6C 

The  Hill  average  of  these  G  =  1/2  (G   +  G  )  is 

^  '  v    r 

1.2[C+C']2  +  5.2CC'  (11) 

4C  +  6C' 

Equation  (11)  is  the  desired  expression  for  the  shear 

constant  in  terns  of  C  and  C1. 

Returning  now  to  the  empirical  relation  of  Fisher  and 

Dever  for  the  shear  constant  C', 

C'  =  0.01391  x  1012  (e*/a)  3*  31t dyne/cm2.        (12) 

This  power  law  relation  for  C1  was  determined  by  taking 
the  best  fit  of  values  of  C1  for  three  Ti-Cr  alloys  and  the 
published  values  of  C'  of  vanadium  and  chromium.   In  this 
fit  the  electron  to  atom  ratio  was  a  parameter.   The  details 
will  be  considered  below.   The  central  result  was  that  by 
assigning  effective  electron  to  atom  ratios,  e*/a,  of  2, 
3,  and  4  to  titanium,  vanadium  and  chromium  respectively 
they  could  obtain  the  relation  (12) . 

Two  important  questions  can  be  raised.   First  why  does 
C1  seem  to  depend  more  on  an  electron  to  atom  ratio,  than 
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on  the  electrostatic  and  ionic  contributions  to  the  shear 
constants  outlined  in  Appendix  II.   Second,  why  the  par- 
ticular numbers  2,3,4  instead  of  the  total  valency  or  some 
other  number.   Quantitative  answers  are  not  forthcoming, 
however  qualitatively  the  reasoning  is  as  follows.   If  in 
the  region  of  (e*/a)  from  2.15  to  4.0,C  is  governed  by  the 
relation  (12)  (the   e*/a  range  of  the  present  work  is  from 
3.0  to  3.46)  then  a  relation  of  the  same  form  and  similar 
magnitude  should  exist  for  the  4d  and  5d  transition  metals, 
Nb,  Ta,  Mo  and  W.   By  plotting  the  published  values  of  C 
C1  and  C  versus  e*/a  for  V,  Nb,  Ta,  Cr,  Mo,  and  W  and  looking 
at  the  percentage  variations  between  the  3d,  4d  and  5d 
elements  they  conclude  that  although  C  does  not  depend 
entirely  on  e*/a  it  has  much  more  of  this  type  dependence 
than  do  Cn  or  C     They  speculate  that  d  electron  ex- 
change interaction  with  nearest  neighbors  could  overcome  the 
negative  contribution  of  ionic  overlap  discussed  in  Appendix 
II. 

For  the  choice  of  e*/a  of  2,3  and  4  for  titanium,  vanadium 
and  chromium  respectively  Fisher  and  Dever  refer  to  a 
recent  APW  band  calculation  of  the  transition  elements  by 
Snow  and  Waber   [Ref .  8] .  In  this  calculation,  Snow  and 
Waber  compute  the  electron  occupancy  of  the  s,  p  and  d 
bands  as  well  as  the  number  of  electrons  in  the  external 
region  between  the  atomic  spheres  that  are  characteristic 
of  the  APW  method.   These  occupancies  are  computed  for  two 
initial  configurations  of  electrons,  s ld      l    and  s2d  where 
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n  varies  from  2  for  titanium  to  9  for  copper.   Snow  and 
Waber  favor  the  S!d  1    configuration  as  the  most  probable 
and,  for  this  configuration,  show  d  shell  occupancies  of 
2.11,  3.07  and  4.01  for  titanium,  vanadium  and  chromium 
respectively.   These  are  very  close  to  n  which  is  the 
e*/a  chosen  by  Fisher  and  Dever  for  the  same  elements.   Thus 
if  the  d  shell  interactions  are  considered  the  primary 
contribution  to  C  this  seems  to  be  reasonably  good  evidence 
for  the  chosen  values  of  e*/a. 

For  the  vanadium-iron  alloy  of  this  work  the  e*/a  of  the 
seven  compositions  could  be  computed  using  the  value  of 
3.0  for  vanadium  and  by  referring  to  Snow  and  Waber ' s  work 
choosing  the  ratio  6.0  for  iron.   They  give  a  d  shell  oc- 
cupancy  of  6.18  for  iron  starting  from  the  s*d    configura- 
tion.  Then  using  equation  (12)  C1  could  be  computed  for 
each  composition.   These  values  of  C  and  the  measured 
values  of  G  can  then  be  used  to  solve  equation  (11) ,  for 
the  shear  constant  C  of  each  composition.   These  results 
are  listed  in  Table  III  along  with  the  best  fit  values  of 
the  derived  quantities. 

In  Figure  6  is  shown  the  results  of  all  the  known  shear 
constant  data  for  the  bcc  region  of  the  first  period  tran- 
sition elements  and  their  alloys.   The  values  of  C  and  C 
shown  for  vanadium-iron  region  are  speculative  in  that  they 
are  derived  from  an  empirical  data  fit  and  a  polycrystalline 
experiment.   However,  the  derived  values  of  C  in  the  reaion 
e*/a  =  3.0  to  3.46  are  continuous  with  the  results  of  the 
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Fisher  and  Dever  alloy  experiment  which  are  in  the  region 
e*/a  =  2.15  to  2.6  of  the  figure.   The  only  other  known 
experiment  in  this  region  is  the  work  on  molybdenum- rhenium 
by  Davidson  and  Brotzen  [Ref.  9]  which  is  shown  in  the  region 
e*/a  =  4.0  to  4.26  of  the  figure.   Fisher  and  Dever  had  no 
explanation  for  how  the  values  of  C  for  vanadium  and  chromium 
joined  up.   These  results  of  the  present  work  merely  state 
that  C  decreases  before  increasing  to  the  chromium  value. 
The  scatter  in  the  C  values  reflects  the  scatter  in  the 
experimental  results.   The  points  for  V-Fe  alloys  on  the  C 
line  in  Figure  6  are  calculated  from  the  empirical  relation 
(12)   and  it  can  be  seen  that  C  apparently  has  a  maximum 
at  e*/a  =  4.0.   Fisher  and  Dever  point  out  that  this  maximum 
corresponds  to  the  maximum  d  band  width  in  Snow  and  Wabers 
calculation. 

Variation  in  the  anisotropy  ratio  with  e*/a  is  readily 
apparent  in  Figure  6.   At  the  left  hand  side  the  ratio  is 
greater  than  unity  and  decreases  to  that  value  for  e*/a 
of  2.8.   It  continues  to  decrease  through  the  V-Fe  region. 
In  the  region  above  e*/a  =  4  the  anisotropy  is  again  in- 
creasing for  the  Mo-Re  alloys.   The  anisotropy  returns  to 
unity  for  e*/a  =  4.26.   Thus  the  concept  of  an  isotropic 
V-Fe  alloy  is  probably  incorrect.   A  more  likely  choice 
would  be  vanadium  with  a  few  atomic  percent  titanium  to 
make  e*/a  =  2.8.   Ti-V  forms  a  bcc  solid  solution  in  this 
region. 
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It  is  interesting  to  note  in  Figure  6  that  C  extrapolates 
to  zero  at  an  e*/a  of  approximately  4.2  which  also  heralds 
the  onset  of  the  sigma  phase  of  the  vanadium-iron  system 
[Ref .  10] .  That  C  would  continue  to  zero  and  create  an  in- 
stability seems  unlikely  and  it  is  thus  assumed  that  C 
must  increase  eventually  toward  the  Cr  value.   In  terms 
of  the  energy  of  an  electron  gas  discussed  in  Appendix  II, 
the  addition  of  electrons  in  the  region  above  e*/a=3.5  would 
push  the  Fermi  surface  closer  to  the  Brillouln  zone  boundary 
where  shear  strain  on  the  direct  lattice,  which  is,  of 
course,  reflected  in  the  reciprocal  lattice  could  increase 
the  average  electron  kinetic  energy  and  contribute  a  positive 
component  to  C.   Another  way  of  visualizing  the  electron 
gas  effect  on  C  is  in  terms  of  the  density  of  states  function 
Snow  and  Waber  show  that  the  density  of  states  at  the  Fermi 
level  varies  with  group  number  in  such  a  way  that  a  minimum 
occurs  for  e*/a  just  less  than  4  for  the  sxd    configuration 
This  is  the  region  where  increasing  e*/a  causes  the  Fermi 
level  to  traverse  between  the  lower  energy  bonding  states 
of  the  d  band  and  the  upper  energy  anti-bonding  states.   The 
correlation  between  these  effects  and  the  requirement  that 
the  C  value  reach  a  minimum  and  turn  up  toward  chromium  is 
unknown.   It  is  only  pointed  out  that  they  both  occur  in 
the  same  region  of  e*/a.   This  concludes  the  discussion  of 
the  shear  constants.   Before  going  on  to  the  bulk  modulus 
a  short  discussion  of  how  the  various  ratios,  e*/a,  were 
determined  is  included. 
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In  making  the  choice  of  e*/a  for  the  elements,  Fisher 
and  Dever  plotted  C'  for  their  three  alloy  compositions, 
vanadium  and  chromium  versus  1,  2,  3  and  4  excess  electrons 
in  chromium  over  those  in  titanium.   Only  the  value  2  gave 
a  smooth  curve  through  the  five  points.   They  then  plotted 
In  C*  versus  In  e*/a  for  and  e*/a  of  Cr  from  1  to  8 .   The 
only  e*/a  to  come  within  the  experimental  error  of  C  for 
their  alloy  and  V  was  e*/a  =  4  for  chromium.   This  fit  was 
1.25%  low  for  the  published  C'  of  Cr.  and  had  the  parameters 
given  in  (12) . 

Since  they  had  established  the  e*/a  for  V  the  remaining 
choice  in  this  work  was  the  value  for  Fe .   Values  of  C  were 
calculated  using  e*/a  =  6  and  7  for  iron.   An  e*/a  of  7 
has  the  effect  of  increasing  the  maximum  e*/a  for  V-Fe  to 
3.61  and  slightly  lowers  the  C  values.   However,  since  G 
is  very  insensitive  to  change  in  e*/a  a  choice  between  6 
and  7  could  not  be  made  on  this  basis  and  it  was  decided  to 
be  consistent  and  use  the  s*d     configuration  with  an  e*/a 
of  6  for  iron.   If  possible  errors  in  the  determination  of 
G  are  considered  it  is  found  that  even  if  G  were  to  increase 
with  alloying  at  the  same  rate  as  the  longitudinal  constant 
the  physical  result  that  C  decreases  with  alloying  for  e*/a 
above  3.0  would  only  be  changed  in  degree. 

Published  values  of  the  bulk  modulus  in  the  transition 
metals  [Ref.  11]  show  that  it  is  peaked  in  each  long  period 
in  the  vicinity  of  Group  VIIIB.   This  dependence  is  similar 
to  that  of  the  cohesive  energy  and  is  independent  of  an 
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interpretation  of  any  model  of  the  contributions  to  the 
bulk  modulus.   Converting  to  electron  concentration,  as  has 
been  discussed  above,  these  data  and  other  published  results 
for  alloys  support  the  idea  that  the  bulk  modulus  increases 
with  increasing  electron  concentration  up  to  a  maximum  at 
e*/a  of  4  or  5. 

One  may  think  of  the  bulk  modulus  as  arising  from  the 
"free"  electron  gas  and  from  the  ion-ion  interactions. 
Snow  and  Waber  using  the  sxd    configuration  estimate  that 
the  s  like  portion  together  with  the  charge  external  to  the 
atomic  sphere  is  between  1.0  and  1.5  electrons.   In  the 
alloys  investigated  here  the  lattice  shrinks  with  the 
addition  of  solute.   This  means  that  if  only  these  electrons 
are  considered  as  contributing  to  the  bulk  modulus  it  will 
increase  with  decreasing  lattice  parameter.   A  simple 
estimate  shows  that  this  contribution  is  about  30%  of  the 
total  for  pure  vanadium  and  roughly  the  same  fraction  for 
the  alloys.   Thus  one  may  infer  that  both  contributions 
increase  with  alloying.   Then  for  the  ionic  contribution 
one  may  conclude  that  alloying  iron  with  vanadium  increases 
the  bulk  modulus  either  because  the  iron-vanadium  ion  inter- 
actions are  larger  or  because  the  iron-iron  interactions 
are  larger.   From  the  shape  of  the  bulk  modulus  curve  the 
latter  interpretation  seems  more  probable. 
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V,   CONCLUSIONS 

Addition  of  iron  to  vanadium  increases  the  bulk  modulus 
but  not  the  shear  modulus,  although  the  latter  data  are 
scattered.   The  increase  in  the  bulk  modulus  can  be  partially 
attributed  to  an  increase  in  electron  gas  concentration  and 
partially  attributed  to  an  increase  in  ion-ion  interactions. 
The  relative  contributions  of  the  s  and  "free"  electrons 
through  the  electron  gas  and  ion-ion  exchange  interaction 
is  approximately  constant  for  the  range  of  alloys  studied. 

Since  the  shear  modulus  remains  constant  this  means  that 
the  shear  constants  C  and  C1  do  not  change  with  alloying  or 
more  probably  that  they  each  change  at  different  rates. 
Using  the  empirical  relation  of  Fisher  and  Dever  for  C 
the  variation  of  both  shear  constants  with  electron  concent- 
ration is  consistent  with  the  previously  measured  values 
for  vanadium  and  the  titanium-chromium  alloys. 

The  anisotropy  ratio  is  a  smoothly  varying  function  over 
the  range  of  effective  electron  concentrations  from  2  to 
4.5.   Thus  this  work  supports  the  idea  of  an  s*d    con- 
figuration for  the  transition  metals  and  an  almost  complete 
dependence  of  C  on  the  electron  to  atom  ratio.   The  concept 
of  an  isotropic  alloy  of  the  transition  elements  has  been 
demonstrated  in  the  Mo-Re  case.   The  existence  of  a  second 
such  alloy  between  two  first  long  period  elements,  V-Ti, 
is  indicated. 
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Future  work  in  this  area  could  include  investigation  of 
the  elastic  constants  of  a  series  of  various  alloys  with 
e*/a  from  2.15  to  4.30.   Such  alloys  as  V-Ti,  V-Cr,  V-Fe 
and  Cr-Fe  might  be  employed  to  investigate  this  entire 
region  of  bcc  stability  in  the  transition  elements. 
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TABLE  I 


No.   Production 
technique 


Vanadium 
source 


Crystal  growth  and 
anneal  technique 


1.    arc  melting 


5 
6 


electron 
beam  melting 


arc  melting 


arc  melting 


electron 
beam  melting 

arc  melting 


arc  melting 


Argonne 

National 

Laboratory 


Bureau  of  Mines 

electrolytically 

refined 

Bureau  of  Mines 
vacuum  arc 
melted 


annealed  I1/2  hours 
at  1000C  longitudinal 
and  transverse  velo- 
city were  measured 
annealed  4  hours  at 
1600C 

float  zone  melting 


strain  anneal 
annealed  2i /z    hours 
at  1600C 

swedged  to  strain 
annealed  I1/2  hours 
at  1600C 


Bureau  of  Mines    passed  slowly  through 
electrolytically   an  induction  coil 


refined 


Naval  Research 
Laboratory 

Argonne 

National 

Laboratory 

Bureau  of  Mines 

electrolytically 

refined 


after  straining 
this  took  4  hours 
at  1600C 

float  zone  melting 
same  as  sample  1 

float  zone  melting 


Notes 

1.  Iron  for  all  samples  was  produced  at  NRL  [Ref.  12] 

2.  Vanadium  produced  by  the  Bureau  of  Mines  is  reported  in 

[Ref.  13  and  14] . 

3.  All  anneals  were  done  in  vacuum 

4.  Temperatures  during  anneals  were  determined  by  optical 
pyrometry . 
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TABLE 

III 

No 

.    Fe% 

P 

3 

G 

e*/a 

C 

C 

B 

A 

V 

0.00 

6.0699 

2.196 

0.474 

3.00 

0.546 

0. 

,431 

1.564 

0.79 

1 

3.25 

6.1809 

2.151 

0.480 

3.10 

0.609 

0. 

,410 

1.511 

0.67 

2 

5.35 

6.2438 

2.230 

0.452 

3.16 

0.649 

0, 

.353 

1.628 

0.54 

5 

7.24 

6.2969 

2.236 

0.490 

3.22 

0.691 

0, 

.388 

1.583 

0.56 

3 

7.40 

6.3017 

2.250 

0.415 

3.22 

0.691 

0, 

.291 

1.696 

0.42 

4 

10.96 

6.3916 

2.284 

0.479 

3.33 

0.773 

0, 

.344 

1.645 

0.45 

6 

13.94 

6.4609 

2.376 

0.483 

3.42 

0.845 

0, 

.326 

1.731 

0.39 

7 

15.32 

6.4917 

2.401 

0.469 

3.46 

0.879 

0. 

.300 

1.776 

0.34 

V 

1 

2 
5 
3 

4 
6 
7 


calculated  best  fit  values 


3.00 

0.546 

0.415 

1.565 

0.76 

3.10 

0.609 

0.388 

1.573 

0.64 

3.16 

0.649 

0.372 

1.588 

0.57 

3.22 

0.691 

0.357 

1.611 

0.52 

3.22 

0.691 

0.357 

1.611 

0.52 

3.33 

0.773 

0.331 

1.672 

0.43 

3.42 

0.845 

0.310 

1.742 

0.37 

3.46 

0.879 

0.301 

1.778 

0.34 

p   in   gm/cm 


[B+^G]  ,    G,    C,    C,    B    in    1012dyne/cm' 
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APPENDIX  I 

The  elastic  theory  of  solids  is  based  on  the  Hooke's 
Law  assumption  that  for  sufficiently  small  strains,  e., 
stress,  T.  is  proportional  to  strain,  where  the  elastic 
constants,  C. .,  are  the  constants  of  proportionality. 

6 

T.=Zc..e-  Al.l 

l   .     11  1 

Deformation  of  the  medium  in  an  adiabatic  manner  results 
in  a  stored  elastic  energy  density,  expressed  by  a  strain- 
energy  function,  w. 

1  6    6 

w  =  -  Z   E   C.e-e.         Al.2 

2  •   •     i~)   i   l 
D=iD=i    J      J 

In  this  notation  the  relation  between  i  or  j  and  the  co- 
ordinate direction  is 

l=xx   2=yy   3  =  zz   4=yz   5=zx   6=xy 
Here  1,  2  and  3  represent  the  normal  strains  which  are 
those  resulting  from  a  stress  on  the  face  of  a  cube  normal 
to  the  x  axis  and  in  the  x  direction  etc.   The  numbers  4, 
5  and  6  denote  the  shear  strains  which  are  those  resulting 
from  a  stress  on  the  face  of  a  cube  normal  to  the  y  axis 
and  in  the  z  direction  etc.   The  shear  strains  in  terms 
of  infinitesimal  displacements  are  given  in  equation 
A1.15.   In  both  cases  the  summations  on  i  and  j  allow  for 
the  general  case  of  an  anisotropic  medium. 
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In  terms  of  the  energy  density  and/or  the  stress  the 
constants  of  proportionality,  C. .,  called  the  elastic 
constants,  stiffness  constants,  or  moduli  of  elasticity 
are 

9T.    .2 

C.  .  =    *   =  I      w    |  (A1.3) 

n   3e .    3e.de.' e=o 

J      D      ID 

It  is  sometimes  convenient  to  express  relationships  in 
terms  of  the  reciprocal  of  the  set  of  equations  (Al.l). 

6 

e.  =  I   S. .  T.  (A1.4) 

1   j=i   ^   D 

The  coefficients  of  this  set  of  equations,  S. .,  are  called 

the  moduli  of  compliance.   Treated  as  matrices  the  product 

matrix  of  the  [S. .]  and  [C. .]  array  is  the  unit  matrix, 

[S.  .]  [C.  .]  =  I. 
ID     ID 

The  strain-energy  function,  w,  is  the  work  done  per  unit 
volume  in  deforming  the  crystal  and  thus  has  the  dimensions 
of  erg-cm3.   The  strains,  e.  are  dimensionless ,  thus  the 

elastic  constants  also  have  the  dimensions  of  an  energy 

—  2 
density  usually  expressed  as  dyne-cm   .   In  this  work  it 


has  been  convenient  to  show  the  C. .  in  units  of  10 12  dyne- 

iD 


cm 

The  number  of  independent  elastic  constants  is  not  36 
as  equation  (A1.2)  might  imply.   From  equation  (Al.3)  and  the 
fact  that  6w  is  a  perfect  differential  we  see  that 

3T.    92w        32w      3T. 

C    =  — i  =  =  =  1  =  c 

ij    3e.    de.de.         de.de.         de .         Si     (Al.5) 
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is  symmetry  C. .=C..  reduces  the  maximum  number  of  indepen- 


Th 

■3    y 

dent  elastic  constants  to  21.   The  elastic  constants  are 
usually  shown  in  a  6  by  6  array  which  for  the  most  general 
case  is 


(A1.6) 


Cll  C12  Ci 3  C14  C15  Ci6 
Cl2  C22  C2  3  Cm  C25  C2G 
Cl  3  C2  3  C3  3  C3^  C35  C36 

Cii»  Cm    Cm    Cm    Cit5  C«f6 

Cl  5  C2  5  C3  5  Cit  5  C55  C5G 
Cl6   C2  6   C3G   CitG   C56   Cg6 

The  number  of  non-zero  independent  elastic  constants  is 
further  reduced,  except  for  the  triclinic  lattice,  by  the 
effect  of  lattice  point  symmetry  operations.   This  effect 
can  be  demonstrated  by  writing  down  the  strain-energy 
function  in  an  original  coordinate  system  and  one  obtained 
from  the  original  by  a  lattice  symmetry  operation.   If  the 
two  strain-energy  functions  are  equated  and  coefficients 
of  like  strains  compared  the  results  are  that  certain 
constants  become  zero  and  others  equal.   For  the  particular 
case  of  cubic  symmetry  only  3  constants  remain  Cn,  Ci2, 
and  Cim.   Thus  the  array  (Al.6)  becomes 
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Cn 

Cl2 

Cl2 

0 

0 

0 

Cl2 

Cn 

Cl2 

0 

0 

0 

Cl2 

Cl2 

Cn 

0 

0 

0 

0 

0 

0 

Cj»  «♦ 

0 

0 

0 

0 

0 

0 

Ci|  if 

0 

0 

0 

0 

0 

0 

c- 

(A1.7) 


Using  this  array  we  can  write  out  equation  (Al.2)  for  the 
strain-energy  function  in  a  cubic  crystal. 


2+e2+e2 

1    2    3 


w  =  —  Ci  i  [e2+e2+e2]+Ci2  [e   e   +e   e   +e   e    ] +-C      [e    2    +e    2    +e      ] 

12  23  312«f«fi+  5  6 


(A1.8) 


A  direct   application   of   equation     (A1.3)   defines  the  cubic 

82w 

elastic  constants  ie  C   =  -s r  etc.   It  is  often  more 

de  i  2 

convenient  when  dealing  with  theories  that  derive  the  elastic 
constants  from  fundamental  concepts  to  define  certain 
combinations  of  the  constants  Cii,  C12  and  Ci»it.   An  equiva- 
lent set  of  these  constants  are 


C  =  1/2  (C11-C12) 

C   =  C«»  it 

B   =  l/3(Cn+  2  C12) 


(A1.9) 


In  this  set  the  constants  C  and  C'  are  shear  constants 

chosen  because  they  can  be  obtained  from  the  strain  energy 

function  by  strains  that  preserve  the  volume  of  the  unit 

cell  to  first  order  in  the  e..   The  constant  B,  is  the 

1 
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bulk  modulus  and  is  the  elastic  constant  corresponding 
to  hydrostatic  compression.   These  three  cases  are  ob- 
tained from  the  cubic  strain-energy  function  (A1.8)  by  choice 
of  the  appropriate  strains. 

1.   C  =  1/2  .(Ci  i-Ci  2 )  •   In  this  case  the  cube  is  compressed 
along  one  axis  and  extended  along  a  second  so  that  two  of 
the  strains,  say,  e   and  e   are  equal  in  magnitude  but 
opposite  in  sign.   All  other  strains  are  zero.   Then 
equation  (Al.8)  becomes 

w  =  (Cii-Ci2)ei2  from  which 


c,=  tt57ler-  (A1-10> 


2.  C  =  C44.   This  corresponds  to  shear  in  a  plane  parallel 
to  a  cubeface  so  that  one  strain,  say  e^  is  finite  and  all 
others  are  zero.   Then  equation  (Al.8)  becomes 

w  =  Ci^eit2  from  which 

92w   1 
C    =  %-%-  (Al.ll) 

3.  B  =  1/3 (C11  +  2  C12)   For  the  case  of  hydrostatic 
compression  the  three  normal  strains  are  equal  and  negative 
and  the  three  shear  strains  are  zero.   Then  equations  (Al.8) 
becomes 

w  =  —  [C11+  2Ci2]ei2  from  which 


2, 
8eV    ei  =  o  (A1.12) 


a  -  —   — — 2 
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A  cubic  crystal  is  said  to  be  isotropic  if  the  two  shear 
constants  in  equation  (A1.9)  are  equal.   The  ratio  of  these 
constants,  A  =  C/C '  is  called  the  anisotropy  and  for  an 
isotropic  crystal  A=l.   Thus,  the  effect  is  to  reduce  the 
number  of  non-zero  independent  elastic  constants  to  two. 
These  are  generally  chosen  to  be  C^  =  y  and  Ci2  -  X. 
The  constants  X   and  y   are  called  the  Lame'  constants. 
The  elastic  constant  Cn  then  maybe  written  in  terms  of 
the  Lame1  constants  and  the  elastic  constant  array  (A1.7) 
contains  only  two  independent  constants. 


(X+2y)         A 

X 

0 

0 

0 

X           (A+2y) 

X 

0 

0 

0 

X                  X 

(X+2y) 

0 

0 

0 

0                  0 

0 

y 

0 

0 

0                   0 

0 

0 

y 

0 

0                  0 

0 

0 

0 

y 

(A1.13) 


Using  equation  (Al.l)  the  isotropic  stress-strain  relations  may 
now  be  written  out. 


Ti  =  (X+2y)e!+  X(e2+e3) 
T2  =  (X+2y)e2+  X(ei+e3) 
T3  =  (X+2y)e3+  X(e2+ei) 
Tu  =  \ieh 
T5  =  ye5 
T6  =  ye6 


XA+2yei 
XA+2ye2 
XA+2ye3 


(A1.14) 
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Here  A  =  ei+e2+e3  is  cabled  the  dilatation.   Following  Mason 
[Ref.  15]  the  expressions  for  the  transmission  of  sound  waves 
in  isotropic  solids  may  now  be  readily  developed. 

In  terms  of  infinitesimal  particle  displacements  u,  v, 
and  w  in  the  x,  y  and  z  directions  respectively  the  strains 
are: 


3u  3v             3w 

ei=  a^  e2=  W      e3=  3^ 

,3v  8w,  ,3u    8w,       ,3v    3u, 

e"=  (yi"  "  37}  es=  (3l  "  ax°  e6=  (a^  ~  ay* 


(A1.15) 


The  equation  of  motion  of  an  infinitesimal  volume  element 
of  density  in  the  x  direction  is 

a2n 
P  hrr   dx  dy  dz  =  F  (A1.16) 

d  t  X 

and  similarly  for  y  and  z.   The  force  Fx  is  that  arising 

from  a  stress  -Ti (x)  on  the  face  at  x  of  a  cube  of  volume 

3T 
AxAyAz  and  a  stress  Ti  (x+Ax)  =  Ti  (x)  +  ■= — Ax  on  the  face 

at  (x+Ax) .   Other  forces  in  the  x  direction  arise  from 

the  shear  stresses  T_  and  T  so  that  in  equation  (Al.16) 

Fv  =  [|£i  +  Sit  +  |Zl]dx  dy  dz  (A1.17) 

x    L3x     3y     3z 

and  similarly  for  F   and  F  . 
J  y       z 

It  one  substitutes  for  the  e.  in  equation  (Al.14) 
equation  (Al.15)  and  take  the  derivatives  indicated  by 
equation  (Al.17)  the  equations  of  motion  (Al.16)  become 
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0=  (x.2y,  fA+p7>u 


32v   .      ,,  .«.;;     9A  • 


tr  =     U+2y)    ^   +    yVzv  (A1.18) 

0  =    (A+2y)    ||  +   yV2w 


2       ,9Z       .3"       .9 


2 


where  v  =    C^r  +  ^pr  +  ^r)     . 

If  the  three  equations  (Al.18)  are  differentiated  with  respect 
to  x,  y  and  z  respectively  and  the  resulting  equations 
summed,  a  wave  equation  for  the  dilatation  is  obtained 


■32A 

9t 


r   =  U+2y)V^A  (A1.19) 


For  a  plane  wave  with  particle  motion  u  only  equation  (Al.19) 
becomes 

P  fp-  =  (A+2y)  |^  (A1.20)* 

which  shows  the  velocity  of  a  longitudinal  wave  traveling 
in  the  x  direction  to  be 


i  p 


vL     ^E 


To  obtain  a  transverse  wave  equation  the  first  and  third 

of  equation  (Al.18)  are  differentiated  with  respect  to  z  and 

x  respectively  and  the  difference  is  formed. 

„  92   ,9u    9w,     n2  r9u    9w,  ,-,  onx 

p  W*    [^  "  9^]  =  ^V   [9z"  -  9^]  (A1'22) 

Then  assuming  a  transverse  plane  wave  traveling  in  the  z 
direction  with  particle  motion  in  the  u  direction, 
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i  (kz-o)t)  /oox  . 

u  =  e         equation  (22)  becomes 


r  =  y  ^rz-r  (Al.23) 


K  3tz    p  9z 
from  which  the  transverse  wave  velocity  is 


S      ||P 


(A1.24) 


Finally,  in  the  discussion  of  polycrystalline  materials 
whose  isotropy  results  from  the  random  orientation  of 
generally  anisotropic  crystallites  the  shear  constant  has 
been  called  G.   Setting  G  equal  to  u  and  using  the  definition 
of  the  bulk  modulus  along  with  the  Lame'  constants 

(X+2y)  =  (b  +  ■*-  G)  (A1.25) 

3 

Substituting  into  (Al.21)  and  (Al.24)  the  polycrystalline  cubic 
material  velocity  equations  become 


VL  = 


[B+^G 

3 


v.-i 


(A1.26) 
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APPENDIX  II 

The  treatment  of  elasticity  in  Appendix  I  is  based  on  a 
classical  Hooke's  Law  force  and  its  attendent  strain-energy 
function.   In  calculating  the  elastic  constants  from  funda- 
mental consideration  one  must  generate  a  function  U  which 
is  the  total  energy  per  atom.   The  equivilance  of  this 
function  and  the  strain-energy  function,  w,  is  assumed  so 
that  in  terms  of  the  energy,  U,  the  definition  of  the 
elastic  constants  is 

C.  .  =  ttM^ — I  (A2.1) 

1 j    de ■ 3e . ' e=o 

The  total  energy  is  a  sum  of  contributions  arising  from 
the  various  forces  between  atoms  in  solids. 

A  successful  theory  of  U  which  is  a  function  of  the 
interatomic  spacing  will  give  first,  the  cohesive  energy 
of  the  crystal  and  secondly  upon  dif f erention,  the  equilibrium 
lattice  constant.   Thirdly,  another  differentiation  of  U 
as  in  (A2.1)  will  give  the  elastic  constants.   This  latter 
requirement  is  a  rigorous  test  of  the  validity  of  any 
fundamental  theory  of  the  total  energy. 

Present  day  theories  of  the  total  energy  U,  which  have 
been  applied  to  the  determination  of  the  elastic  shear 
constants  C  and  C,  originated  with  Fuchs .   Fuchs  applied 
the  methods  of  Wigner  and  Seitz  for  determining  the  cohesive 
energy  of  metalic  sodium.   He  extended  the  treatment  to 
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copper,  and  later  calculated  the  elastic  constants  of 
copper  [Ref .  1] . 

Fuch ' s  calculation  was  applied  to  transition  elements 
by  Zener  [Ref.  16]  and  Isenberg  [Ref.  17].   They  concluded 
that  the  observed  stability  of  the  bcc  transition  elements 
resulted  largely  from  a  repulsive  interaction  between  second 
nearest  neighbors  arranged  in  a  anti-ferromagnetic  array. 
Neighbors  and  Smith  use  Fuch ' s  calculation  to  discuss  dilute 
alloying  of  6ther  elements  with  copper  [Ref.  18]. 

It  is  instructive  to  consider  the  calculation  of  Wigner 
and  Seitz  for  the  cohesive  energy  of  metalic  sodium  since 
the  contributions  to  the  total  energy  for  this  case  can  be 
extended  to  other  monovalent  metals. 

The  electronic  configuration  of  sodium  is  a  single  s 

electron  outside  of  a  closed  neon  core.   In  the  solid 

state  this  electron  is  assumed  to  be  free  to  move  through 

the  crystal  lattice  and  Wigner  and  Seitz  calculate  the  energy 

of  the  lowest  state,  U  ,  of  the  free  electron.   The  cal- 

o 

culation  is  facilitated  by  imagining  the  lattice  to  be 
divided  into  polyhedra  constructed  by  erecting  planes  which 
are  perpendicular  bisectors  of  the  lines  joining  an  atom 
to  its  neighbors.   Such  polyhedra  are  adjacent  and  fill 
all  space  occupied  by  the  lattice.   The  motion  of  the 
electron  through  the  lattice  then  becomes  a  motion  between 
adjacent  identical  polyhedra.   Such  a  motion  will  be  period- 
ic with  respect  to  the  dimensions  of  the  polyhedra  and  the 
boundary  conditions  on  the  electron  wave  function  will  be 
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that  the  function  be  continuous  between  polyhedra  and  that 
it  approach  the  zone  boundary  normally.   On  the  average 
a  monovalent  metal  will  have  only  one  electron  per  polyhedron. 
This  is  a  result  of  both  Coulomb  repulsion  and  the  effects 
of  the  Pauli  exclusion  principle.   Thus  each  electron  finds 
itself  at  a  given  instant  confined  to  a  atomic  polyhedron 
influenced  by  the  potential  of  the  central  ion  core.   Wigner 
and  Seitz  assumed  a  potential  by  Prokofjew  [Ref.  19]  to 
calculate  the  wave  function  and  U  ,  the  lowest  state  energy 
of  this  one  electron  system. 

Since  the  atomic  polyhedron  is  not  very  different  from 
a  sphere,  they  assumed  that  it  may  be  replaced  by  a  sphere 
of  equal  volume.  The  radius  of  this  sphere  is  r  ,  defined 
by  —  iTr  3=Q,   where  Q   is  the  atomic  volume. 

*   3     O 

The  radial  wave  function  R  =  r^ (r)  so  that  the  bound- 
's \y 

ary  condition  on  ¥  which  is  t-      =  0  becomes 
2  9r'r=r 

o 

9Ri      =R  (A2.2) 

8r'  r=r   r 
o 

A  convenient  way  to  obtain  the  dependence  of  the  energy 
on  the  radius  of  the  sphere  is  to  solve  the  radial  Schrodinger 
equation  for  the  wave  function  at  arbitrary  energies  U  . 
After  obtaining  the  wave  function,  R(r) ,  a  tangent  drawn 
from  the  origin  to  the  curve  will  satisfy  the  boundary 
condition  (A2.2)  at  two  radii.   Thus  from  each  radial  function 
two  points  on  the  U(r)  curve  are  determined.   Calculation 
of  several  radial  function  allows  the  U(r)  curve  to  be 
mapped  out.   For  sodium  such  a  curve  has  a  minimum  and  the 
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lattice  would  be  stable  without  considering  other  contribu- 
tions to  the  energy. 

In  fact  however,  the  electrons  in  the  metal  are  dis- 
tributed over  states  above  U  ,  the  energy  of  the  lowest 
state,  by  the  Fermi  distribution  function.   Thus  to  the 
energy  of  the  lowest  state  must  be  added  the  average  ad- 
ditional electron  kinetic  energy.   This  energy  is  deter- 
mined from  the  fermi  energy 

where  the  average  kinetic  energy  of  an  electron  is 

UF  =  j  Up  (A2.4) 

These  are  the  primary  contributions  to  the  cohesive  energy 
of  sodium.   However,  sodium  has  a  relatively  small  ion 
core  size  so  that  the  increase  in  the  energy  of  the  lattice 
from  core  overlap  forces  is  not  significant.   The  larger  ion 
core  size  in  copper  does  not  allow  this  assumption  and  the 
repulsive  exchange  energy  between  ions  must  also  be  con- 
sidered.  This  energy  is  usually  assumed  to  have  the  form 

Ux   =  |^  (A2.5) 

where  W  is  the  repulsive  bond  energy  of  a  single  pair  of 
ions  and  S  is  the  lattice  coordination  number.   W  is  usually 
taken  to  be  of  the  form 

W  =  De"r/p  (A2.6) 
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Since  p  is  a  small  fraction  of  the  interatomic  distance 
only  first  and  second  nearest  neighbors  are  usually  con- 
sidered to  take  part  in  the  interaction.   Other  energies 
which  are  small  and  do  not  contribute  substantially  to 
the  cohesive  energy  are  the  Van  der  Waals  energy  and  a 
correction  to  the  lowest  state  energy.   This  correction 
UT  is  the  result  of  polyhedron-polyhedron  interaction  and 

Li 

is  the  sum  of  two  terms. 

UL  =  UE  -  Us  (A2.7) 

where  \J      is  the  energy  of  a  lattice  of  positive  point 
charges  immersed  in  a  negative  electron  sea  and  U   in  the 
energy  of  a  positive  point  charge  in  a  sphere  of  negative 
charge. 

The  fact  that  a  certain  energy  contributes  signifi- 
cantly to  the  cohesive  energy  does  not  imply  that  it  will 
significantly  affect  the  elastic  constants.   In  obtaining 
the  elastic  shear  constants  of  copper  from  the  energies 
outlined  above  Fuch ' s  used  only  the  coulomb  attraction 
energy  of  U„  in  equation  (A2.7)  and  the  overlap  repulsion 
energy  of  equations  (A2.5)  and  (A2.6). 

The  reasons  the  other  energies  do  not  contribute  to 
the  shear  constants  are  as  follows: 

a)   For  the  equilibrium  lattice  spacing  U  (r)  is  close 
to  an  inflection  point  and  its  second  derivative 
is  close  to  zero. 
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b)  For  monovalent  metals  the  Fermi  surface  is  near 
the  center  of  the  Brillouin  zone  and  is  approxi- 
mately spherical.   Thus  shear  strains  do  not  alter 
the  average  electron  energy. 

c)  The  Fermi  energy  and  the  spherical  contribution 
to  U  ,  and  Uc ,  are  volume  dependent  energies. 

Li  O 

However,  the  strains  chosen  to  define  C  and  C* 
conserve  the  volume  and  thus  these  energies  do 
not  contribute  to  the  to  the  elastic  constants. 

d)  The  Van  der  Waals  energy  makes  only  a  minor 
contribution  to  C  and  C1. 

Thus  Fuch's  retained  only  two  contributions  to  the  energy 

U  =  U   +  IT  (A2.8) 

1     hi 

where  U^,  is  the  electrostatic  potential  energy  of  a  lattice 
of  point  changes  +  e  imbedded  in  a  uniform  negative  change 
distribution  of  -e  per  atomic  volume  and  the  energy  U   is 
the  exchange  repulsion  energy  between  closed  shell  ion 
cores . 

For  a  bcc  lattice  UF  can  be  evaluated  using  the  Ewald 
method  [Ref.  20]  to  do  the  lattice  sum  and  the  electro- 
static contributions  to  the  elastic  constants  are  then 
calculated  using  equation  (A2.1)  so  that 

e2  e2 

C'   =  0.0997  —       C^,  =  0.7423  -7        (A2 . 9 
E  a1*        E  a* 

where  a  is  the  lattice  constant. 
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Similarily  ionic  stiffnesses  C'    and  C   are  given  by  Fuch ' s 
in  terms  of  the  energy  W  of  equation  (A2.6) 

CI   =   -N    [rW*]       +   ^N    [r2W  ']    +  ±®    [rW] 

I  3  112  z  2  2  22 

CT    =   -N    [r2W'']    +   -5-N    [rW']       +   N    [rW']  (A2.10) 

1  9  119  11  2      2 

Here  allowance  has  been  made  that  the  nearest  neighbor  in- 
teraction W  may  differ  from  the  next  nearest  neighbor 

interaction  W    Primes  indicate  differentiation  and  the 
2 

subscripts  on  the  brackets  denote  evaluation  at  a  distance 
r  between  ions  corresponding  to  first  or  second  nearest 
neighbors . 

The  total  shear  constants  are  then 


C  =  CT  +  C_       C*  =  Cl  +  C'  (A2.ll) 

In  practice  C   and  Cl  are  usually  determined  by  cal- 
culating C„  and  C'  and  subtracting  this  from  the  measured 
constants.   Featherston  and  Neighbours  [Ref.  21]  have  done 
this  for  V,  Nbf  Ta,  Mo,  W  and  Fe .   They  note,  in  agreement 
with  Isenberg,  that  for  all  cases  except  Fe  the  term 

—  N(rw')    in  Cl   is  negative  and  so  must  be  offset  by  the 
3      1  1       I 

second  nearest  neighbors. 

The  application  of  this  analysis  to  dilute  alloys  of 
copper  was  done  by  Neighbors  and  Smith.   By  introducing 
parameters  to  accomodate  the  effects  of  solute  atoms  on 
the  electron  to  atom  ratio  and  the  ionic  overlap  contribu- 
tion and  by  correcting  for  changes  in  the  stiffnesses  due 


54 


to  changes  in  the  lattice  constant  with  solute  they  calculate 
corrections  to  C  and  C. 

Further  extensions  of  Fuch ' s  theory  to  the  transition 
elements  and  their  alloys  has  encountered  the  difficulty 
of  dealing  with  the  status  of  the  d-electrons .   Such 
electrons  are  more  tightly  bound  than  s  electrons,  they 
thus  occupy  a  relatively  narrow  band  with  a  correspondingly 
higher  density  of  states  than  s  electrons.   Further,  the 
d  wave  functions  are  more  directional  in  character.   The 
multi-valency  of  these  metals  and  the  changes  in  e/a  with 
alloying  will  certainly  cause  changes  in  the  fermi  surface 
which  could  result  in  an  electron  gas  contribution  to  the 
shear  constants.   The  antif erromagnetic  array  that  Zener 
and  Isenberg  predicted  using  Fuch ' s  work  has  only  been 
observed  experimentally  for  chromium. 
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